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(N ■ Abstract. A generalized variant of the Calderon problem from electrical 

(-H ' impedance tomography with partial data for anisotropic Lipschitz conduc- 

tivities is considered in an arbitrary space dimension n > 2. The follow- 
ing two results are shown: (i) The selfadjoint Dirichlet operator associated 
with an elliptic differential expression on a bounded Lipschitz domain is deter- 
mined uniquely up to unitary equivalence by the knowledge of the Dirichlet- 
to-Neumann map on an open subset of the boundary, and (ii) the Dirichlet 
operator can be reconstructed from the residuals of the Dirichlet-to-Neumann 
map on this subset. 
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1. Introduction and main results 

Let C be a uniformly elliptic and formally symmetric second order differential 
expression of the form 

<N _ _ 

(1.1) C = - d jajk d k + J2 ( a i d 3 ~ d ^j) + a 

j,k=l j=l 

■ with variable coefficients on a bounded Lipschitz domain f2 in R™, n > 2. The main 

objective of the present paper is to show that the selfadjoint Dirichlet operator 

(1.2) Au = jCu, domA= {u e H^(fl) : Cue L 2 (n)} , 

associated with C in L 2 (f2) is uniquely determined by a local variant of the Dirichlet- 
to-Neumann map on some open subset ui of the boundary d£l, and that A can be 
reconstructed from the residuals of this partial Dirichlet-to-Neumann map. Here 
the Dirichlet-to-Neumann map M(X) is defined by 

(1.3) M(A) : H'/^dn) -> H-'/^dn), u x \ 9n ^ d c u x \ dn , 

where u\ 6 H (0) is a solution of Cu = Xu and dcu\\gn denotes the conormal 
derivative of u at dfl. The mapping M(X) is well-defined for all A in the resolvent 
set p(A) of A; see Section 2 for further details. 

The above inverse problem is closely connected to and inspired by the famous 
Calderon problem from electrical impedance tomography, where the aim is to de- 
termine the isotropic conductivity 7 of an inhomogeneous body from current and 
voltage data measured on the surface or on parts of it. The classical Calderon 
problem corresponds to the special case ajk = 7#jfc, a,j ■ = a = 0, 1 < j, k < n, in 
(1.1), and the knowledge of M(A) on d£l or cj C d£l for some A or at A = 0; see 
[20]. 

The Calderon problem has been a major challenge in the field of inverse prob- 
lems for PDEs in the last three decades. The first positive results were obtained for 
uj = dQ using only Dirichlet and Neumann data for A = in the pioneering work 
[67] in dimension n > 3 and smooth 7, see also [59] for 7 G C 1,1 (fi) and [56] for 
the reconstruction of 7 from the boundary measurements. In the two-dimensional 
case the first main contribution was the solution of the problem for 7 £ W 2,p {Vl) 
in [57]; later in [10] conductivities in L°°(0) were allowed. For partial data given 
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only on special subsets of dfl uniqueness was shown in the recent works [19, 45] 
for a C 2 -function 7 and a reconstruction method was provided in [58], see also [41] 
for a generalization in the two-dimensional case. Also the more general case of 
an anisotropic conductivity (ffljk)"fc = i has been investigated; in this situation, the 
single coefficients are in general not uniquely determined. Nevertheless, unique- 
ness up to diffeomorphisms was first shown for real-analytic coefficients assuming 
knowledge of M(0) in [51] in dimension n > 3 and in [66, 68] for n = 2; more 
general cases were treated in [9, 27, 65]. In the publications [49, 50] the related 
problem of determining a real-analytic Riemannian manifold from the given Dirich- 
let and Neumann boundary data on arbitrary open subsets of dQ was considered. 
There uniqueness up to isometry in n > 3 and uniqueness of the conformal class 
in n = 2 was shown for partial data supported in an open subset of d£l; see also 
[13, 14, 15, 27, 43, 44, 47] for closely related problems as, e.g., the multidimensional 
Gclfand inverse spectral problem and inverse problems for the wave equation with 
elliptic data. For a detailed recent review and further references we also refer to 
[69]. 

The aim of the present paper is to prove somewhat different, milder types of 
uniqueness and reconstruction results in space dimension n > 2 for partial data 
given on an open subset uj of dQ. Since the coefficients of C are not uniquely 
determined in general we focus on the selfadjoint Dirichlet operator (1.2) associated 
with C on f2. In return this point of view onto the problem allows us to consider 
the more general differential expression (1.1) and to impose, in an arbitrary space 
dimension n > 2, the following comparatively weak assumptions on the coefficients 
of C. 

Assumption 1.1. The coefficients ajk and aj are bounded Lipschitz functions on 
fl satisfying cikj = ajk, 1 < j,k < n, and a G L°°(f2) is real- valued. Moreover, C is 
uniformly elliptic, i.e., 

n n 

E >cE** 
j,k=i fe=i 

holds for some C > 0, all £ = (£1, . . . £„) T S M™ and x G O. 

Furthermore, we impose the following conditions on the domain f2 and the subset 
to of the boundary <9£1 where Dirichlet and Neumann data is assumed to be given. 

Assumption 1.2. il is a bounded Lipschitz domain in K™, n > 2, and w is a 
nonempty, open subset of the boundary dSl. 

The following two theorems on uniqueness and reconstruction of the Dirichlet 
operator A in (1.2) from the knowledge of the Dirichlet-to-Neumann map onwC dfl 
are the main results in the present paper. The usual duality between H 1 / 2 (di}) and 
H^^idn) is denoted by (•,•)• 

Theorem 1.3. Let ft anduj be as in Assumption 1.2, and let £1 and £ 2 be differen- 
tial expressions on Q, of the form (1.1) with coefficients djk,i, aj,i, ai andajk,2, a j,2, 
a 2 satisfying Assumption 1.1. Denote by A\, A 2 , and Mi, M 2 the corresponding 
Dirichlet operators and Dirichlet-to-Neumann maps, respectively. Assume that 

(Mi(\)ip,ip) = (M 2 (AV, V) fie ff 1/2 (dft), suppi^, ip C uj, 

holds for all A G T>, where T> C p(Ai) fl p(A 2 ) is a set of points which has an 
accumulation point in p{A\) n p(A 2 ). Then A\ and A 2 are unitarily equivalent. 

We point out that Theorem 1.3 is a "mild" uniqueness result in the sense that 
even in the special case ajk — 7^fe, fflj = a = 0, 1 < j, k < n and oj = dQ 
in the classical Caldcron problem it does not imply uniqueness of 7 as shown in 
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[10, 19, 45, 56, 57, 58, 59, 67]. Under additional smoothness assumptions Theo- 
rem 1.3 can also be viewed as a consequence of the gauge equivalence of second 
order elliptic operators on manifolds determined from their boundary spectral data 
in the multidimensional Gelfand inverse spectral problem; cf. [13, 15, 43, 44, 47]. 

In order to state our second main result we remark that A M> M(A) is a meromor- 
phic operator function; the residual of the Dirichlet-to- Neumann map M at a pole 
/i is a linear operator from 7? 1 / 2 (9i7) to H~ x / 2 (dVl) which is denoted by Res^M. 
Since the aim in the next theorem is to use only knowledge of M and its residuals 
on u!, also the operator Res M M is replaced by a local version. The residual Res^M 

on lo is defined on the linear subspace (d£l) of functions ip <G H 1 ^ 2 (dn) with 
supp (p C to by 

Res^M : H^dSl) -»■ (l^/ 2 (df>))', (Res^M^ty) := (Res M M^>), 

1 /2 

where ip, ip £ HJ (<9f2). The next theorem provides a one-to-one correspondence 
between the eigenfunctions of the Dirichlet operator A and the residuals of the 
Dirichlet-to-Neumann map M on w and, hence, yields a representation of A in 
terms of these residuals. 

Theorem 1.4. Let Assumption 1.1 and Assumption 1.2 be satisfied and let M(A) 
be the Dirichlet-to-Neumann map associated with the elliptic differential expression 
C. Then the poles of A i-> M(A) coincide with the eigenvalues of the Dirichlet 
operator A and for each eigenvalue Afe of A the corresponding eigenspace is given 
by 

kcr(A - A fe ) = ran (r^ 1 o Res" fc M) , 

where Tfc denotes the restriction of the Neumann trace operator u M> dcu\ u onto 
kcr(A — Afe). In particular, there exist ip^\ . . . , <p^) k \ £ H 1 ^ 2 (dfl) with support in 
w, such that the functions 

ef } := Tl : l {Re^M)pf\ i = 1, . . . , n(k), 

form an orthonormal basis ofkei(A — Afe) and A can be represented in the form 

oo n(k) 

Au = Afe ( u > e i fe) )i,2(fi) 4 k \ u G domA. 

k=l i=l 

The proofs of our uniqueness and reconstruction results are based on the powerful 
interplay of modern operator theory with classical PDE techniques, as, e.g., unique 
continuation results from [7, 8, 40, 70]. Some of the main ideas are inspired by 
abstract methods in extension and spectral theory of symmetric and selfadjoint 
operators as provided in [2, 11, 12, 22, 23, 24, 25, 42, 48]. Further operator theoretic 
approaches to elliptic boundary value problems and related questions via Dirichlet- 
to-Neumann maps or other analytic operator functions can also be found in the 
recent publications [3, 4, 5, 16, 17, 18, 31, 33, 34, 35, 36, 63, 64]. For general methods 
from extension theory of symmetric operators that are applied to elliptic PDEs we 
also refer to, e.g., [6, 30, 32, 37, 53, 54, 60, 61], the monographs [29, 38, 52, 55], 
and the references therein. 

The main part of the present paper is devoted to the proofs of the two main 
results; along the proofs we show Proposition 2.5 and Proposition 2.7 on the eigen- 
values and eigenspaces of the Dirichlet operator, as well as a density result in 
Lemma 2.6 which is of independent interest. The paper closes with a short appen- 
dix, which summarizes some basic facts on unbounded operators in Hilbert spaces 
and on Banach space- valued analytic functions. 
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2. Proofs of Theorem 1.3 and Theorem 1.4 

In this section we give complete proofs of the uniqueness and reconstruction 
theorems from above. Instead of two single proofs the material is ordered in several 
smaller statements which then lead to the proofs of the main results. 

We fix some notation first. By H s (ft) and H s (dtt) we denote the Sobolev spaces 
of order s > on Q and dSl, respectively, and by Hq(Q) the closure of the set of 
C°°-functions with compact support in il with respect to the iP-norm. Further, 
H~ s (d£l) denotes the dual space of H s (d£l); the duality is expressed via 

(2.1) /(?) = (/, 50 = fo7), f&H~ s (dct), ^eip(do), 

which extends the L 2 inner product. For the Lipschitz domain we write u\an G 
H^idd) for the trace of u g H 1 ^) at the boundary dfl and d c u\ dn g H- 1 / 2 {dQ) 
for the conormal derivative or Neumann trace of it at dVl (with respect to the 
differential expression £), see, e.g., [55, Chapter 4] for more details. Recall further 
that Hq(ZI) coincides with the kernel of the trace operator u t-> u\dn on _ff 1 (17). 
In order to define a restriction of the conormal derivative to the nonempty, open 
subset uj C dCl let 

(2.2) H^ 2 (dQ) = {<pe ff 1/2 (<9ft) : supp if C uj} 

be the linear subspace of H 1 / 2 (dfl) which consists of functions with support in ui. 
The restriction dcu\ u of the Neumann trace dcu\on to uj is defined as 

(2.3) (d c u\ u ){<p) := (d c u\ m )(<p) = (d c u\gn,<p), V G H^\dQ). 

Let us recall some well-known properties of the Dirichlet operator associated to 
C, which can be found in, e.g., [29, Chapter VI] and [55, Chapter 4]. 

Proposition 2.1. The Dirichlet operator A in (1.2) is a selfadjoint operator in 
L 2 (il) and its spectrum cr{A) consists of isolated (real) eigenvalues with finite- 
dimensional eigenspaces. The Dirichlet eigenvalues accumulate to +oo and are 
bounded from below. 

The next lemma shows that the Dirichlet-to-Neumann map and the Poisson 
operator in Definition 2.3 below are well-defined. 

Lemma 2.2. For all A in the resolvent set p{A) of A and all tp g H 1 / 2 (dQ) there 
exists a unique solution u\ g H l {£l) of the boundary value problem 

(2.4) Cu = \u, u\an = f- 

In particular, for all A g p(A) the Dirichlet-to-Neumann map M(A) in (1.3) is 
well-defined. 

Proof. Let A g p(A) and <p g H 1 / 2 (dfl). Then the homogeneous problem 

(£ - X)u = 0, u\gn = 0, 

has only the trivial solution, and by [55, Theorem 4.10] it follows that the inhomo- 
geneous problem (2.4) has a unique solution. □ 

Besides the Dirichlet-to-Neumann map a Poisson operator which maps functions 
on dfl onto the corresponding solutions of (2.4) will play an important role. 

Definition 2.3. Let A g p(A). The Poisson operator is defined by 

(2.5) 7 (A) : H^ 2 (dn) -> L 2 (n), u x \ an -> u x , 

where u\ G i? 1 (fi) is the unique solution of (2.4) with ip = u\\qq. The range of 

1 /2 

the restriction of the Poisson operator to Hj (dCl) is denoted by j\f\, 

(2.6) Af\ = {u g iJ^O) : Cu = Au, supp(u| a n) C uj} . 
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The operator 7(A) is well-defined for each A G p(A) by Lemma 2.2 and the 
relation M(X)ip = dc{'j{X)y>)\dQ. holds for all <p G H 1 / 2 (dfl). Some properties and 
formulas for the Poisson operator and the Dirichlet-to-Neumann map will be given 
in the next lemma. Its proof is essentially based on the second Green identity 

(2.7) (Cu, v) - (u, Cv) = (u\ d n, d c v\ d n) - (d c u\ d n, v\an) 

for u,v G if 1 (£2) satisfying Cu,Cv G L 2 (£l), see, e.g., [55]. Here (•,•) on the left 
hand side denotes the inner product in L 2 (il) and on the right hand side the duality 
between H^ 1 ^ 2 (dfl) and H 1 ^ 2 (dil); cf. (2.1). In the following it will be clear from 
the context whether the entries in (•, •) are functions on f2 or dSl, respectively, so 
that no confusion can arise. We remark that in a more abstract setting statements 
of similar form as in Lemma 2.4 can be found in, e.g., [11, Proposition 2.6], [23, 
§ 1], and [48, § 2]. 

Lemma 2.4. Let A, fx G p{A), let 7(A), 7(/z) be the Poisson operators and let M(A), 
M{p) be the Dirichlet-to-Neumann maps. Then the following statements (i)-(iv) 
hold. 

(i) 7(A) is bounded and the adjoint operator 7(A)' : L 2 (Q) — > H^ 1 / 2 (dQ) is 
given by 

7(A)'u - -d c ((A - A)-^) |an, u G L 2 (n). 

(ii) 7(A) and 7(/x) satisfy the identity 

7(A) = (7+(A- m )(A-A)- 1 ) 7 (m)- 

(hi) The Poisson operators and the Dirichlet-to-Neumann maps are connected 
via 

(77 - A) ( 7 (AV, 7^) = (M(X)ip, V) - {<p, M(M), 

and, in particular, (M(\)(p,ip) = {ip,M{\)ijj) holds for alltp,ip G H 1 / 2 (dfl). 
(iv) M(A) is bounded, the function A i-» M(A) is holomorphic on p{A), and the 
identity 

(M(AV,V) = (^,M(A )V) 

+ (A - A) ((J + (A - A )(A - A)- 1 ) 7 (Ao)(^,7(Ao)^) 

/10/ds /or a/? A, Ao G p(^4) fl^rf V, ^ G H 1 ' 2 (dil). In particular, every 
eigenvalue of A is either a pole of order one or a removable singularity of 
the mapping A i-> M(A). 

Proof, (i) In order to verify the formula for 7(A)' we show first 

(j(\)ip,u) — (<p, —dc((A — ~\)~ 1 u)\dii) for all ip G H 1 ^ 2 (dfl), ueL 2 {Q). 
Let p G ff 1 / 2 ^), set u A = 7 (A)v? and let u G L 2 {9.). One finds by (2.7) 
( 7 (AV, «) = (u A , (/ + \(A - A)- 1 ) u) - (Au A , (A - A)" 1 **) 
= (u A , £(A - A) _1 «) - (r« A , (A - A) _1 «) 
= {dcu\\dn, (A - A) _1 M|ao) - (uA|an,3c ((A - A) _1 u) | an ) 
= (<p, -dc ((A - A) _1 u) |an) , 

where we have used (A — A) _1 u|an = in the last equality. Then it follows with 
(2.1) from 

(7(A)'u)(¥>) = W,7(AV) = («,7(A)v) = - (dc({A - A) _1 m) |an) (<p) 
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that 7(A)' acts as in the assertion and is defined on L 2 (Q). Moreover, the above 
reasoning also implies that 7(A) is closed and hence bounded by the closed graph 
theorem. 

(ii) For X,p G p(A), if e i? 1/2 (<9Sl) and u G L 2 {tt) wc find by (i) 

( 7 (AV, u) - ( 7 ( M V, u) = (p, -d c ((A - m)" 1 ^ - Jl)(A - A)- 1 n)|ao) 
= ( 1 (p)<p,(X-p)(A-X)- 1 u) 

= ((\-v)(A-\y 1 1 (^,u). 

Hence we have j(X)ip = 7 (ju)v + (A — p)(A — A) _1 7 (/x)y>, which shows (ii). 

(iii) Let X, p G p{A) and ip,ip G H x / 2 {dVt), and set = 7(A)^s and u M = j(p)tp. 
Again by (2.7) we find 

(7*-A)(7(AV,70u)V) - {u x ,£v„) - - (M(A)<^) - (p,M(/i)^). 

(iv) From (M(A)v?, ^) = {<p, M{X)ip) for A G p{A) and (p,ip <E H^ 2 {dn) it follows 
that M(A) : H 1 / 2 ^) ->• H- 1 ^ 2 {dn) is closed and hence bounded by the closed 
graph theorem. Furthermore, by (iii) we have 

(M(A)p, V) = (<f, M(Ao)V) + (Ao - A) ( 7 (A)<^, 7 (A )^) , 

which together with (ii) implies the formula in (iv). The remaining statement in 
(iv) follows from this and the corresponding properties of the resolvent of A, see 
also Appendix. □ 

The statement in Lemma 2.4 (iv) on the singularities of M and the eigenvalues 
of A will be improved later in Corollary 2.8, where it turns out that the Dirichlet 
eigenvalues coincide with the poles of the meromorphic operator function A h-» M(A) 
and its restriction on cu. 

The following proposition is essential for the proofs of our main results. It states 
that the restriction of the Dirichlet operator A onto {u G dom A : dcu\u = 0} is an 
operator without eigenvalues. In an operator theoretic language this implies that 
this restriction is a simple symmetric operator. The main ingredient in the proof 
of Proposition 2.5 is a classical unique continuation theorem for solutions of second 
order elliptic differential inequalities, see, e.g., [70] and [7, 8, 40]. Moreover, the 
proof makes use of the sesquilinear form <&/;(•, •) induced by C on _ff 1 (fi), 

r, n n 

&c(u, v) := I f 2, a-jkdku ■ djV + \ , ( a j(dj u ) ' v + ' djv) + auvjdx 
■>'• ^ 1 i=i 1 

for u,v G _ff 1 (0), and the corresponding first Green identity 

(2.8) $c(m, v) = (Cu, v) + (d c u\ dn , v\ an ) 

for all u, v G iJ 1 (SI) satisfying Cu G L 2 (VL)\ see [55, Chapter 4]. 

Proposition 2.5. There exists no eigenfunction u of the Dirichlet operator A 
satisfying dcu\ u — 0. 

Proof. Let SI D SI be a bounded Lipschitz domain such that <90 \ oj C 9S1 and 
contains an open ball O. We extend the coefficients ajk, dj, and a of the differential 
expression C from (1.1) to functions a,jk, o,j, and a on SI such that Assumption 1.1 
holds for the corresponding differential expression C in SI defined as in (1.1). 

Assume that there exists A and u 7^ in the domain of the Dirichlet operator (1.2) 
with Cu — Xu on O and d^u]^ — 0. Since we have u\qq — 0, we can extend u by 
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zero on £1 \ Q to a function u £ Moreover, Cu £ L 2 (Q) and Cu = Xu holds. 

In fact, we compute for tp £ C^°(f2) 

(Cu - Xu) (tp) = (u, <p ) — X(u, tp ) , 

where the left hand side is understood in the sense of distributions and the right 
hand side consists of integrals on f2 with integrands vanishing outside of f2. De- 
noting the restriction of <p to Q by tp, it follows with the help of the first Green 
identity (2.8) that 

(Cu - Xu)(tp) = §c(u,w) - X(u,lp) = {Cu,Jp) + (d c u\ dn ,Tp\ dn ) - X(u,lp). 

Since Cu — Xu, supp(^|ao) C w and dcu\ w = we conclude together with (2.3) 

(Cu- Xu)(tp) = (d c u\ dn ,7p\ 9n ) = (dcu\ w )(tp\ d n) = 0. 

Hence Cu £ L 2 (VL) and Cu = Xu hold; in particular, u is locally in H 2 , see, e.g., [55, 
Theorem 4.16]. Furthermore, we obtain 

n n n n 

- ^2 a-jkdjdku = (^X - a + ^ d 3 a^ju + ^ ^ ^ dja jk -a k + a k ^j (d k u). 

j,k=l j=l k=l j = l 

Since the functions a,j k and cij, 1 < j, k < n, together with their derivatives of first 
order as well as a are bounded on il and u = on fl \ there exist constants a 
and (3 such that 



^ cij k djd k u 



<a\u\+pY,\ d ^\ 

k=l 



(2.9) 

3,k=l 

holds a.e. on O. As u = on O it follows from the differential inequality (2.9) and 
classical unique continuation results that u vanishes identically on O; cf. [70]. In 
particular, we conclude u = on ft, a contradiction, since u was chosen to be an 
eigenfunction of A. □ 

Our last preparatory lemma will establish, as a consequence of Proposition 2.5, a 
density statement on the ranges Af\ of the Poisson operators 7(A) in (2.5) restricted 

1 /2 

to Hj (dfl). Recall that Af\ is the space of solutions u of the boundary value 
problem (2.4) which satisfy supp(«|an) C w. 

Lemma 2.6. Let O C C+ be an open set and let O* = {X £ C : A G O} . Then 

spanjA^A : A e OUO*} = span{ 7 (A)(^ : tp £ H^ 2 (dn), X £ OUO*} 
is a dense subspace of L 2 (il). 

Proof. The proof consists of three separate steps. It makes use of two further 
operator realizations S and T of the differential expression C. We consider the 
restriction 

Su = Cu, domS = {u £ Hq(CI) : Cu £ L 2 (Q), d c u\ u = 0}, 

of the Dirichlet operator A in L 2 (il), which has no eigenvalues by Proposition 2.5, 
and we define the operator T in L 2 (f2) by 

Tu = Cu, domT = {u £ H\n) : Cu £ L 2 {Q), u\ dn £ H^ 2 (dQ)}. 

It follows immediately that the Dirichlet operator A is a restriction of T and that 
the spaces A/a coincide with ker(T — A). In the first step of this proof we show that 
these spaces are dense in the spaces ker(5 f * — A). In the second step, which can also 
be found in a different form in [46], a selfadjoint restriction of S in the orthogonal 
complement of span {A/a : A £ O U O* } is constructed. In the last step we then show 
that the spectrum of this selfadjoint operator is empty, which implies the assertion. 
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Step 1. In this step we show that Af\ = kcr(T — A) arc dense subspaces of 
ker(S* - A), A G O U O* . For this we check first that S = T* holds. In fact, for 
u G domS* and v G domT the second Green identity (2.7) together with u\en = 0, 
dcu\u = 0, and supp(u|on) C u> implies 

(Tv, u) = (v, Cu) + (v\dn, dcu\on) ~ (dcv\on,u\dn) = (v, Cu); 

cf. (2.3). Hence u G domT* and T*u = Cu = Su by the definition of the adjoint 
operator. For the converse inclusion let u G domT*. From A C T we obtain 
T* C A* = A and therefore T*u = Cu and u G domA In particular, we have 
u G H x (£i), Cu G L 2 (fl), and u\on = 0. It remains to show dcu\ w = 0. For 
v G domT we have supp(w|gn) C to and from (2.3) and (2.7) we obtain 

(d c u\ u )(v\ dn ) = (d c u\ d n,v\ d n) = -(T*u,v) - (u,Tv) + (u\ dn , d c v\ d n) = 0. 

1 /2 

As runs through HJ (dft) as v runs through domT, it follows dcu\ u = 0, 
hence u G dom5*. We have shown T* = S. This implies T = T** = S* and hence 
the spaces N\ — ker(T — A) are dense in the spaces ker(S'* — A). 
Step 2. In this step we show that the Hilbert space 

(2.10) M := (span{AA A : A G O U 0*})^ 

is invariant for S and that Sm S \ (A4 n domS 1 ) is a selfadjoint operator in A4; 
cf. [46]. Observe first that by step 1 the symmetric operator £ = T* is closed and 
hence ran (S — A) is closed for all A G C \ K; cf. Appendix. According to step 1 we 
also have 

M = (span{kcr(5** - A) : A G OUO*}) 1 = p| ran (S -A). 

xeouc 

Let u G M n domS". Then u G ran (5" - A) for all A G O U O* , i.e., for each 
A G O U O* there exists u\ G domS 1 such that (S — X)u\ = u holds. This implies 

Su = S(S - X)u x - (S - X)Su x G ran (S - A) 

for each A G O U £>* and hence Su G M. Therefore, S M = S \ (M n domS 1 ) is 
an operator in and since S is symmetric we conclude that the restriction Sm is 
a symmetric operator in M. For the selfadjointness of Sm m A4 it is sufficient to 
show 

(2.11) ran (Sm — v) = A4 = ran (5Vi — v) for some v G C 

For this let u G M. We claim that w := (S — ^) _1 u G domS 1 belongs to M. In 
fact, for each A G O U O* , A 7^ u, the function 

«A ^^((S-A)- 1 -^-^)- 1 ),, 

satisfies (S — X)v\ = v and hence v G ran (S — A) for all A G U O* , A ^ In order 
to check that also v G ran (5 — v) holds, we choose a sequence (A n ) n G 0, \ n ^ t/, 
which converges to v. As above it follows that (S — \ n )~ 1 u G ran (S — A) for all 
A G O U 0*, A ^ A„, and, in particular, (£ — \ n )~ 1 u G ran (S — v). Since 5* is a 
closed symmetric operator the estimates ||(5 — ^) _1 || < |Imi^| _1 and ||(5 — A„) _1 || < 
ImA„| _1 hold, and hence 

v-(S- Xny'u =(S- v)-^ - (S - A„)" 1 « ={v- X n )(S - v)- 1 ^ - X n )~ 1 u 
implies 

(2.12) v = lim (S - A„) _1 m. 

n— >oo 

Since (S— \ n )~ 1 u G ran (S — v) and ran (5* — ^) is closed we conclude v G ran (S — v) 
from (2.12). We have shown v G M. Aloreover, (Sm — v)v = (S — v)v = u and, 
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hence, the first equality in (2.11) holds. The second equality in (2.11) follows in 
the same way. Therefore Sm = S \ {M fl dom5) is a selfadjoint operator in M. 

Step 3. It follows from step 2 that the operator S can be written as the direct 
orthogonal sum Sm © S M ± with respect to the decomposition L 2 (Q) = A4 (B A4 J ~, 
where Sm is a selfadjoint operator in M. and S M ± = S \ (A4 1 - ndom S) is a closed 
symmetric operator in M.^~. Moreover, the selfadjoint Dirichlet operator A admits 
the decomposition A = Sm ©^i, where A' ± is a selfadjoint extension of Sm 1 - 
Since the spectrum of A consists only of eigenvalues (see Proposition 2.1) the same 
holds for the spectrum of the selfadjoint part Sm- Clearly, each eigenfunction of 
Sm is also an eigenfunction of £ = Sm (£> Sm 1 - i but by Proposition 2.5 this operator 
has no eigenfunctions. Therefore, the spectrum of the selfadjoint operator Sm in 
M is empty which implies M — {0}; cf. Appendix. Hence 

L 2 (fl) = M 1 - = span{A/A : A e OllO*} 

by (2.10), where span denotes the closed linear span. This completes the proof of 
Lemma 2.6. □ 

With these preparations we are ready to prove the uniqueness result Theorem 1.3. 

Proof of Theorem 1.3. Let C\ 1 C 2 be elliptic differential expressions as in the 
theorem and let Ai, A 2 and Mi, M 2 be the corresponding selfadjoint Dirichlet 
operators and Dirichlet-to-Neumann maps, respectively. The associated Poisson 
operators from Definition 2.3 will be denoted by 71(A) and 72(A), respectively, 
A G p{Ai) n p{A 2 ). Since V C p(Ai) n p(A 2 ) in Theorem 1.3 is a set with an 
accumulation point in the intersection of the domains of holomorphy of the functions 
Mi and M 2 (see Lemma 2.4 (iv)), it follows that Mi and M 2 coincide on oj, i.e., 

(2.13) (Mi(AV,V) = (M 2 (AV,V), tp^eH^idQ), 

holds for all A e p(Ax) n p(A 2 ), and, in particular, for all A e C\M. Here Hl /2 (d£l) 
is the linear subspace of iJ 1//2 (9fi) which consists of functions with support in cj; 
cf. (2.2). 

Let in the following X, p E C\R, \^=~p- Lemma 2.4 (iii) and (2.13) yield 

/ m m a (Mi(A)y,^)-(y,Mi(^) 

(2 U) fi 

= = — t = (72(A) v, 72^)^) 

fl — A 

for all ip, ip e Hl ,2 (dVL). Next we define a linear mapping V in L 2 (f2) on the domain 
domV = span{7i(AV : V G #i /2 (<9ft), A e C \ K} by 

t t 

V : ^7i(A J -)^ ^ ^72(A j )^-, A, 6C\l, ft € H^/ 2 {dil), 1 < j < 1 

Formula (2.14) yields that V is a well-defined, isometric operator in L 2 (il) with 

1/2 

ran V = span{7 2 (A)<p : (p <E (90), A 6 C \ R}. By Lemma 2.6 the domain and 
range of V are both dense subspaces of L 2 (tt), and hence the closure U of V in 
L 2 (il) is a unitary operator in L 2 (f2). Obviously {/7i(A)<p = 7 2 (A)<£ holds for all 

1 /2 

A e C \ M and p e (dCl). With the help of Lemma 2.4 (ii) one computes for 
A^M 

U( Al - A)-S^ - U ^ X \- U ^ = (A 2 - A)- 1 ^^ 

A — jl 
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1 /2 

for tp £ Hj (dQ,). From Lemma 2.6 we then conclude 

UiAi-X)- 1 = (A 2 -\)- 1 U 

and therefore UAiu — A 2 Uu for all u G domAi, that is, A\ and A 2 are unitarily 
equivalent. □ 

The next proposition is the key ingredient in the proof of the reconstruction 
result Theorem 1.4. Here the residual Res A M : H^ 2 (dn) -> H^^idVt) of the 
Dirichlet-to-Neumann map at a pole A and a local version of it play the main role. 
Define the operator 

(2.15) Res^M : i^ /2 (<90) -> (i^ 2 ^))', 95 i-> Res^M^, 

where Res^-M^ is the restriction of ResA-Mtp G H^ 1 ^ 2 (dfl) onto H^ 2 (dQ,), i.e., 

(2.16) (Rcs^M^)W) := (Res A M^)(V>) - (Res A M^), i> G i^ /2 (dft); 

cf. the definition above Theorem 1.4 in the introduction. Note that in the special 
case lo = dfl the operators Res^M and Res^M coincide. It turns out that the 
Neumann trace operator on uj maps each eigenspace of the Dirichlet operator bi- 
jectively onto the range of the corresponding residual of the Dirichlet-to-Neumann 
map on ui. 

Proposition 2.7. For each A the mapping 

t\ : ker(j4 — A) — > ranRes^M, u i-> dcu\ u , 
is an isomorphism and, in particular, 

ran Res^M = {tp:th ere existsu G ker(A — A) such that dcu\ u = ?/>}. 

As an immediate consequence of this proposition we obtain the following state- 
ment which complements Lemma 2.4 (iv). 

Corollary 2.8. A point X is an eigenvalue of A if and only if A is a pole of 
the Dirichlet-to-Neumann map on uj. In this case the dimension of the eigenspace 
kcr(A — Ao) coincides with the dimension of the range of the operator Rcs A M. 

Proof of Proposition 2. 7. Note first that if A is not an eigenvalue of A, then the 
function M is holomorphic at A by Lemma 2.4 (iv) and hence the residual ResA-M 
is zero. Therefore the assertion holds in this case. 

In the following let A (g R) be an eigenvalue of A and let T\u — dcu\ u be the 
Neumann trace operator on u on ker(A — A) defined in (2.3). For u G kcr(^4 — A) 
and some fixed /ieC\lwe obtain for all ip G H^ 2 (dfl) by Lemma 2.4 (i) 

(tau)(V0 = (dculenrf) = (d c {{A-ji)- l Au- (A-Jiy^u) | an ,^) 

<2 ' ' " - ((A - ~p)dc ((A-Jtr'u) |an,^) = (CP - A) 7 (m)V?)- 

Let us show that t\ is injective. Assume there exists u G ker(A— A), u 7^ 0, such 
that T\u = 0, that is, 

(2.18) (d c ({A--p)- 1 u)\^) = Q 

for all V G Hl ,2 {dVL) by (2.17). Since u ^ 0, also (A - ^ holds, but 

(A- \)(A-Jl)- 1 u= (A-JI)-\A- X)u = 

implies that (A — ~p)~ 1 u is an eigenfunction of A with vanishing Neumann trace on 
uj by (2.18). This contradicts Proposition 2.5. Hence t\ is injective. 

Next we show that t\ maps onto ranRes^M, which by (2.17) is equivalent to 

(2.19) ranRes"M = ran f ker(yl - A)). 
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The inclusion C in (2.19) will be shown first. For this denote by E\ the spectral 
projection onto the eigenspace ker(A — A) of A. Let r\ G C \ R, r\ ^ Jl, and O be an 
open ball centered in A such that O does not contain any further eigenvalues of A 
and n, T) £ O. Denote by Y\ the boundary of O and write the spectral projection 
E\ as a Cauchy integral; cf. Appendix. For if, ip e H 1 / 2 {dti) we obtain with the 
help of Lemma 2.4 (ii) and (iii) 



2ttz 
1 
2?ri 



— / (J- 



-(7(0^,7^)^) 



iff 



1 



C-v 



(MCO^VO t n 

(^-CXm-C) i'P-Oi'P-v) {(-v)(~p-v) 



(M(J7)^,V) 



cf. [26, §1.1]. Note that the second and third fraction in the integral are holomorphic 
in O as functions of ( and hence can be neglected. In the remaining fraction, we 
develop the function ( n- M(£) into a Laurent series at A. Since this function has 
a pole of, at most, order one in A, see Lemma 2.4 (iv) and the Appendix, we obtain 

= ^/r, (C-A) (?7 -C)(M-C) dC = (,-A)(7l-A)' 

where Cauchy's integral formula was used in the second equality. Therefore we 
have 

(Res\Mcp, ip) = (v~ A) (71 - A) (E xl ( V )^ j(^) 
and hence, in particular, 

(2.20) Res^M^ = (7 7 -A)(7Z-A) 7 (M) / ^A7(^, ^ G l^ /2 (dO); 

cf. (2.15) and (2.16). This shows the inclusion C in (2.19). 

In order to show the inclusion D in (2.19) let u 6 kcr(A — A) and e > 0. Since 
■y(fj,)'E\ is bounded by Lemma 2.4 (i), there exists 5 > such that \\u — v\\ < S 
implies \\j(fj,)'E\u - j(fi)'E\v\\ < e. Note that E\u — u as u 6 kcr(A - A). By 

Lemma 2.6 we find £ G N, r]j G C \ R, and ^ G Hl /2 (dfL), l<j< £, such that 



i=i 



< 5. 



Then we have 
(2.21) 



< e 



and from (2.20) and (2.21) we conclude 



7(a*)'"-5Z 



Res^M^j 



ivj -A)(m-A) 



< e, 



that is, 7(/i)'w € ranRes^M. The fact that ker(A — A) is finite-dimensional 
(see Proposition 2.1) shows (2.19). Thus t\ is bijective and maps the finite- 
dimensional space ker(A — A) onto ranRes^M, i.e., t\ is an isomorphism. □ 



Theorem 1.4 is now essentially a consequence of Proposition 2.7. 



Proof of Theorem 1.4. For the statement on the poles of M see Corollary 2.8. 
In order to prove the representation of A let us denote by (Xk)k the sequence of 



12 



JUSSI BEHRNDT AND JONATHAN ROHLEDER 



eigenvalues of A and by Tfc the restriction of the Neumann trace operator onto the 
eigenspace ker(A — Afe) from Proposition 2.7, 

r k = r Xk : kcr(A - Afe) -> ranRes^ fe M, u i-)- <9£u| w . 

Since by Proposition 2.7 Tfc is an isomorphism, the formula 

kcr(A - A fe ) = ran (r^T 1 o Res^ fe M) 

for the eigenspace corresponding to the eigenvalue Afe follows immediately. In par- 
ticular, we can choose ip[ k \ . . . , y&k) e Hl/ 2 (d£l), n(k) = dimkcr(^4 — Afe), such 
that the functions 

ef^r-^Res^M)^, 1 < i < n(k), 

form an orthonormal basis in ker(A — Afe). Then the orthogonal projection E k in 
L 2 (Q) onto ker(A — Afe) is given by 

n(fe) 

Since the spectrum of A consists only of eigenvalues with finite-dimensional eigen- 
spaces we conclude that A can be represented in the form 

oo n(k) 

Au = J2 X k^2(u,e < i k) )e < i k) , u E dom A 

fe=l i=l 

□ 

Remark 2.9. Since M is a holomorphic operator function whose values are bounded 
operators from if 1//2 (<9S!) — > H^ 1 ^ 2 (dQ) it is sufficient to assume knowledge of 
M (X)ip in Theorem 1.3 and Theorem 1.4 for ip in an arbitrary dense subspace of 

1 /2 

ffi' (<9Q) and A in a discrete set of points T> with an accumulation point in p(A). 
In particular, since the spectrum of A is discrete, it is possible to choose V as a 
sequence in any nonempty, open subset of R. 

Appendix 

In this appendix we recall some definitions and basic facts on (unbounded) op- 
erators in Hilbert spaces and Banach space-valued mappings, which can be found 
in, e.g., the monographs [1, 21, 28, 39, 62]. 

Linear operators in Hilbert spaces. Let T-L and /C be Hilbert spaces with scalar 
products (•, -)u and (•, -)k, and let T be a linear operator from H to JC. We write 
domT, kerT and ranT for the domain, kernel and range of T, respectively. The 
operator T is said to be closed if its graph is a closed subspace of H x K,. If domT 
is dense in T-L, then the adjoint operator T* is defined by T*g :— g, where 

domT* = {g E K, : exists g e H such that {Tf,g) K = {f,g) H for all / E domT}. 

Observe that T* is a closed linear operator from JC to %. Moreover, domT* is 
dense in tC if and only if the closure T of (the graph of) T is an operator, and in 
this case T** coincides with T. 

Let A be a linear operator in T-L. Then A is said to be symmetric if the relation 
{Af,g)u = (f,Ag)-u holds for all f,g E dom A If A is densely defined, then A 
is symmetric if and only if the adjoint operator A* is an extension of A, that is, 
dom .4 C dom A* and A* f = Af holds for all / e dom A; in short A C A*. If 
A = A* holds, then the operator A is called selfadjoint. Recall that a symmetric 
operator A is selfadjoint if and only if ran (A — A±) = H holds for some X± E C ± 
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and that for a closed symmetric operator A and A e C\l the operator (^4 — A) -1 
defined on the closed subspace ran(A — A) is bounded by |ImA| _1 . 

Let S be a closed linear operator in H. The resolvent set p(S) consists of all 
A G C such that (S — A) -1 is a bounded operator defined on %. Note that the 
set p(S) is open in C. The spectrum <r(S) of S is the complement of p(S) in C, 
in particular, o~(S) contains the eigenvalues of S, i.e., the points A £ C such that 
ker(S' — A) ^ {0} holds. Recall that for a symmetric operator A in H the eigenvalues 
are real and that for a selfadjoint operator A the whole spectrum a(A) is contained 
in R. Moreover, by the spectral theorem the spectrum of a selfadjoint operator A 
is empty if and only if the Hilbert space T-L is trivial. 

Holomorphic functions with values in Banach spaces. Let X be a Banach 
space and let D C C be an open set. If the function m : D — > X is holomorphic on 
D and p £ C is a pole of m, then the residual of m at p is 



where L M denotes a closed Jordan curve with interior O which contains p, such that 
C \ {lA Q D. Equivalently, Res M m is the first coefficient of negative order in the 
Laurent series expansion of m at p. If A is a selfadjoint operator in a Hilbert space 
H, then A H> -Ra(A) = (^4 — A) -1 is a holomorphic function on p(A) with values in 
the space of bounded linear operators in %. Here each isolated point p in <r(A) is a 
pole of Ra of order one and the orthogonal projection E u in H onto the eigenspace 
kei(A — p) is 
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